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We consider the motion of a particle on a surface which is a small perturbation of the standard 
sphere. One may qualitatively describe the motion by means of a processing great circle of the 
sphere. The observation is employed to derive a subsidiary Hamiltonian system that has the form of 
equations for the top with a 4-th order Hamiltonian, and provides the detailed asymptotic description 
of the particle’s motion in terms of graphs on the standard sphere. 
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I. INTRODUCTION 


The dynamics of a particle which is allowed to move freely, i.e. without the action of external forces, on a smooth 
surface is the classical problem in analytical dynamics, Q. It is generally hard to solve. In fact, the four Hamiltonian 
equations for orbits on a surface can be reduced to a system of two Hamiltonian equations by use of the integral of 
energy and elimination of time, 0. But the system obtained in this way may have no further integrals and admit of no 
exact solutions. Thus, the usual reduction method, Q , or the momentum map according to the current terminology, 
does not work. The problem looks even more pessimistic if one aims at drawing a picture of the ensemble of orbits 
on a surface, for it requires the study of the general form and disposition of orbits on a surface of general shape. 

In this paper the surface is supposed to be a perturbed standard sphere. Using the perturbation theory we construct 
a Hamiltonian system, which enables us to give a fairly detailed picture of the ensemble of orbits by means of graphs 
on the standard sphere; the vertices of the graphs corresponding to orbits which are asymptotically closed and the 
edges of the graphs to orbits joining the almost closed ones. 

To be specific, the central idea relies on the circumstance that if the surface does not differ substantially from a 
sphere, great circles of the latter may serve a good approximation to orbits of the particle. If segments of an orbit 
are short enough, one can visualize it as winding up in coils; loops or rings of the coil corresponding to great circles 
of the sphere, see FIGID Hence approximating the successive rings by great circles, we may describe the change in 
the position of the rings by the precession of a great circle, which in its turn is determined by the normal vector L 
of the plane cutting the sphere by the great circle. To cast this picture in a more quantitative form we may use the 
fact that the normal vector L is the angular momentum of the particle moving on the great circle. Thus, we may use 
perturbation theory, i.e. averaging method, for determining the motion of the angular momentum. 


II. AVERAGED EQUATIONS OF MOTION 


The equations of motion of a particle on a surface given by the equation ip{x) = 0 can be cast in the form of the 
equation. 
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FIG. 1: Two coils I, II of an orbit on the surface ipix) = 0; vectors Ni and Nu are the normals to the planes of great circles 
approximating the coils 


X = X 


dip 

dx 


( 1 ) 


The Lagrangian multiplier can be found explicitly , so that the equation of motion, in the form that does not involve 
A, reads 


1 , ^ 

dx 

\dx) 

In this paper we assume that the equations of the surface be of the form 


( 2 ) 


^{x) = '^{xi +eixf) - I, 
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where £i are small. Then equations © read 
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The angular momentum L = f x p. verifies the equations, which follow from 0 






3 


Li 

1/2 


-4 


-4 


-4 


Ej(2 + 12£jx])x/ 

Ey(2xj +4ejx]y 

Ej(2 + 12£jX^)xE 
Ej(2a:j +4ejx3)2 

Ej (2 + 12 £ j - a :: 2)^-^.2 

Ej(2a;i +4eja;3)2 


a;2a;3(e3a;3 


X3Xi(£iXi 


XiX2(£2X2 


£2X2) 


£3X3) 


£1X3). 


( 4 ) 


The equations given above are exact, in the sense that they do not involve any approximation and do not use the £i 
being small. Their treatment still needs further refining, and this will be done with the help of the method of averaging. 
Generally, the approach relies on studying the evolution equations for integrals of motion of the unperturbed system, 
i.e. in our case the normals to the planes of the large circles, with respect to the basic periodic solution of the latter. 
The averaging serves as a filter separating the main regular part of the solution from the oscillating one caused by 
small terms considered as perturbation, see |^. 

We may write the basic equation for the particle’s motion on the sphere of unit radius in the form 


X = cos{ujt + 9)ei + sin{ujt + 6 )e 2 


where vectors ei, 62 , £3 can be conveniently taken in the form 
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ei = 


62 = 


\/L2 + 
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L^/L. 
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{-Ll- LlLiL2,LiL3) 


£3 — J-{Li,L2, L3). 


. ‘2 

The angular velocity lo is given by the equation uj'^ = x = valid to within the first order of perturbation. Here 
Li,L 2 ,L 3 are coordinates of the normal to the plane of the great circle determining the solution, i.e. the angular 
momentum. 

Let us turn to the exact equations for the angular momentum 0- With the help of the equations given above and 
neglecting terms of the second, and higher, order in the £i, we can transform equations 0 in the form 
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It should be noted that the right-hand sides of equations comprise terms oscillating in time and terms that 
vary slowly. By using the averaging method, 0, that is on neglecting the oscillatory terms we obtain the averaged 
equations for the angular momentum 


Ll = - —[(£3 - e2)Li -I-£3L2 - £2L3] , 

L2 = ^ +{ei - e 3 )L\+eiL"^ , ( 6 ) 

L3 = ^ ^2^ [^^LJ - £iL 2 -I-(£2 - £i)L 3 ] . 


It is worth noting that equations have the Hamiltonian form determined by the usual Poisson brackets for the 
angular momentum, |^, 

{Lj, Lj'\ — ^ ^ ^ijk^k^ 

k 

and the Hamiltonian 



(7) 


This circumstance is particularly interesting because, usually, the averaging procedure is not compatible with Hamilto¬ 
nian structure. The system we have obtained is the integrable Hamiltonian one, but its exact solution is cumbersome. 
Therefore, we shall find a qualitative description of the system’s motion and extensively use numerical simulation. 

The important point is considering the stationary solutions to equations for which the right-hand sides turn out 
to be zero. They correspond to trajectories that are closed to within oscillating terms neglected during the averaging. 
Their equations split into three parts SI, S2 and S3, determined by conditions on Si, as follows. 
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FIG. 2: Separatrix net corresponding to the graph of Type I on the sphere. 


S3 Vector L subject to Lio ^ 0, L 20 7 ^ 0, L 30 7 ^ 0 and the quadratic constraints imposed on ec 

r2 r2 r2 

_ -^10 _ _ _ -^^20 _ _ _ ^30 _ 


£1 £2 ~ £2 £3 + £3 £1 £1 £2 + £2 £3 ~ £3 £1 ~£i £2 + £2 £3 + £3 £1 


It is worth noting that equations S2 involve the fulfilment of the inequalities e 2£3 > 0, £361 > 0, and £i £2 > 0 for 
cases S2.a, S2.b, S2.c, respectively, whereas equations S3 involve 


£1 £2 - £2 £3 + £3 £1 > 0 
£1 £2 + £2 £3 - £3 £1 > 0 
—£1 £2 + £2 £3 + £3 £1 > 0 


Linearizing equations m at the stationary solutions and, considering small fluctuations of L round them, we may 
study their stability, which turns out to be determined by the requirements 

51 a. £i£2 > 0; 
b. £2£3 > 0 ; 

C. £3£i > 0. 

52 a. £i£2 - £2£3 + £3£1 < 0; 
b. £i£2 + £2£3 - £3£i < 0 ; 

C. -£i£2 + £2£3 + £3£i < 0. 

53 any Si . 

We may put these equations in a more graphic form by using the integral = const, and consider the motion 
of Z on a sphere of fixed radius, the integral of energy H taking appropriate values. Then the stable solutions are 
fixed points as regards equations ®, the stable and the unstable points are foci and saddle points, respectively, the 
separaterixes being lines joining the fixed points. Together, they generate a graph on the sphere, having the fixed 
points as vertices and the separaterixes as edges. It is important that the separaterixes, i.e. the edges of the graph, 
are oriented according to time, t, so that the graph is the oriented one, and invariant with respect to the symmetry 
R —R and time inversion t —t. 
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FIG. 3: Type I trajectories of the auxiliary system on the semi-sphere; dashed lines are typical trajectories, the solid ones 
separatrixes. Parameters of deformation: e\ = 0.02, £2 = 0.03, £3 = 0.04. 



FIG. 4: Type II trajectories of the auxiliary system on the semi-sphere; dashed lines are typical trajectories, the solid ones 
separatrixes. Parameters of deformation: £1 = 0.01, £2 = 0.03, £3 = 0.04. 


Now we are in a position to determine the graphs by employing the computer simulation of the equations ®, 
and using the types of the fixed points found above. It should be noted that we must check as to whether the 
solutions provided by equations m agree with those given by original equations m , see FIGli The phase picture 
can be obtained by constructing a mesh generated by solutions to equations taking into account the types of fixed 
points. The results are illustrated in FIGsEHi 

We obtain the following types of the graphs: 


Type I FIG El 7 foci and 6 saddles; being subject to the constraints: 

£i£ 2 - £2£3 + SsSl > 0 

£i£2 + £2£3 - £3£i > 0 . ( 8 ) 

~£ i £2 + £223 + £ 3 £i > 0 

Type II FIG0 5 foci and 4 saddle points; £i are not equal to zero, have the same sign, and at least one of equations 
is not true. 
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FIG. 5: Type III trajectories of the auxiliary system on the semi-sphere; dashed lines are typical trajectories, the solid ones 
separatrixes. Parameters of deformation: ei = —0.02, £2 = 0.03, £3 = 0.04. 



FIG. 6 : Type IV trajectories of the auxiliary system on the semi-sphere; dashed lines are typical trajectories, the solid ones 
separatrixes. Parameters of deformation: £1 = —0.01, £2 = 0.00, £3 = 0.01. 


Type III FIGEI 3 foci and 2 saddle points; Si being subject to one of the following constraints: £2^3 > 0 and ei£2 < 0 ; 
£3£i > 0 and e2£3 < 0; £162 > 0 and £3£i < 0. 

Type IV FIGEI 2 foci and 1 saddle point; Ei being subject to one of the following constraints: £1 = 0 and £2£3 < 0; 
£2 = 0 and £ 3 £i <0; £3 = 0 and £i£2 < 0. 


The separatrix nets depend on values of the coefficients of the deformation Si, and generate regions I, II, III, IV in 
the £i space. 

It is important that the lines dividing the domains corresponding to types I and II, FIG[ 7 | are given by the 
homogeneous equations 


1 . 

2 . 

3 . 


£i£2 - £2£3 + £3£i = 0 
£ i £2 +£2£3 - £3^1 = 0 
— £ i£2 + £2£3 + £3^1 = 0 


( 9 ) 
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FIG. 7: Regions of Si corresponding to Types I - IV of the solutions to the auxiliary system. The white area indicates Type 
I solutions, the filled and the barred ones Type II and Type III. The lines dividing the Type I and Type II regions are subject to 
equations 




FIG. 8: Comparison of solution to: A. the initial equations for geodesics; B. the averaged equation given by the auxiliary 
system. 


The type of a graph corresponding to the solutions is completely determined by the numbers of foci and saddle 
points. The dependence of the conformations of the foci and the saddles on values of Si is illustrated in FIGEI 
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III. CONCLUSION 


The main instrument of the present investigation is the auxiliary Hamiltonian system, which can be considered 
as a reduction of the initial problem to a dynamical problem on specific configuration and phase spaces. Points of 
the new configuration space are geometrical objects, i.e great circles, of the configuration space, i.e. the standard 
sphere. Thus, we obtain a Hamiltonian system that describes the transformation of these objects. In a sense, the 
approach follows the classical method of Klein and Lie, 0 , of constructing a new space with objects of the given 
one. For the specihc case of an ellipsoid close to the standard sphere our classification of orbits are in agreement with 
the classical results by Jacobi, for geodesics on ellipsoid. We feel that our asymptotic approach is useful for the 
treatment of more general problems, for example the motion of rigid bodies, and intend to consider the problem in 
the subsequent paper. In analytical terms we study an asymptotic reduction of the system of equations for orbits 
on a deformed sphere to that of the top, but with the Hamiltonian of the fourth order. The simplification we get in 
this way, is substantial. Indeed, the initial Hamiltonian system could be non-integrable,whereas the auxiliary system 
is totally integrable and described by a graph that comprises vertices, which correspond to stationary solutions, or 
almost closed orbits, and edges, which can be visualized as orbits joining them, that is orbits which continuously 
approach more and more closely to coincidence with the closed ones. 
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